Predicting lncRNA-protein interactions through deep learning framework employing multiple features and random forest algorithm.

RNA-protein interaction (RPI) is crucial to the life processes of diverse organisms. Various researchers have identified RPI through long-term and high-cost biological experiments. Although numerous machine learning and deep learning-based methods for predicting RPI currently exist, their robustness and generalizability have significant room for improvement. This study proposes LPI-MFF, an RPI prediction model based on multi-source information fusion, to address these issues. The LPI-MFF employed protein-protein interactions features, sequence features, secondary structure features, and physical and chemical properties as the information sources with the corresponding coding scheme, followed by the random forest algorithm for feature screening. Finally, all information was combined and a classification method based on convolutional neural networks is used. The experimental results of fivefold cross-validation demonstrated that the accuracy of LPI-MFF on RPI1807 and NPInter was 97.60% and 97.67%, respectively. In addition, the accuracy rate on the independent test set RPI1168 was 84.9%, and the accuracy rate on the Mus musculus dataset was 90.91%. Accordingly, LPI-MFF demonstrated greater robustness and generalization than other prevalent RPI prediction methods.

Deep Conservation and Unexpected Evolutionary History of Neighboring lncRNAs MALAT1 and NEAT1.

Long non-coding RNAs (lncRNAs) have begun to receive overdue attention for their regulatory roles in gene expression and other cellular processes. Although most lncRNAs are lowly expressed and tissue-specific, notable exceptions include MALAT1 and its genomic neighbor NEAT1, two highly and ubiquitously expressed oncogenes with roles in transcriptional regulation and RNA splicing. Previous studies have suggested that NEAT1 is found only in mammals, while MALAT1 is present in all gnathostomes (jawed vertebrates) except birds. Here we show that these assertions are incomplete, likely due to the challenges associated with properly identifying these two lncRNAs. Using phylogenetic analysis and structure-aware annotation of publicly available genomic and RNA-seq coverage data, we show that NEAT1 is a common feature of tetrapod genomes except birds and squamates. Conversely, we identify MALAT1 in representative species of all major gnathostome clades, including birds. Our in-depth examination of MALAT1, NEAT1, and their genomic context in a wide range of vertebrate species allows us to reconstruct the series of events that led to the formation of the locus containing these genes in taxa from cartilaginous fish to mammals. This evolutionary history includes the independent loss of NEAT1 in birds and squamates, since NEAT1 is found in the closest living relatives of both clades (crocodilians and tuataras, respectively). These data clarify the origins and relationships of MALAT1 and NEAT1 and highlight an opportunity to study the change and continuity in lncRNA structure and function over deep evolutionary time.

Structural basis for inactivation of PRC2 by G-quadruplex RNA.

Polycomb repressive complex 2 (PRC2) silences genes through trimethylation of histone H3K27. PRC2 associates with numerous precursor messenger RNAs (pre-mRNAs) and long noncoding RNAs (lncRNAs) with a binding preference for G-quadruplex RNA. In this work, we present a 3.3-Å-resolution cryo-electron microscopy structure of PRC2 bound to a G-quadruplex RNA. Notably, RNA mediates the dimerization of PRC2 by binding both protomers and inducing a protein interface composed of two copies of the catalytic subunit EZH2, thereby blocking nucleosome DNA interaction and histone H3 tail accessibility. Furthermore, an RNA-binding loop of EZH2 facilitates the handoff between RNA and DNA, another activity implicated in PRC2 regulation by RNA. We identified a gain-of-function mutation in this loop that activates PRC2 in zebrafish. Our results reveal mechanisms for RNA-mediated regulation of a chromatin-modifying enzyme.

Hit identification of novel small molecules interfering with MALAT1 triplex by a structure-based virtual screening.

Nowadays, RNA is an attractive target for the design of new small molecules with different pharmacological activities. Among several RNA molecules, long noncoding RNAs (lncRNAs) are extensively reported to be involved in cancer pathogenesis. In particular, the overexpression of lncRNA metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) plays an important role in the development of multiple myeloma (MM). Starting from the crystallographic structure of the triple-helical stability element at the 3'-end of MALAT1, we performed a structure-based virtual screening of a large commercial database, previously filtered according to the drug-like properties. After a thermodynamic analysis, we selected five compounds for the in vitro assays. Compound M5, characterized by a diazaindene scaffold, emerged as the most promising molecule enabling the destabilization of the MALAT1 triplex structure and antiproliferative activity on in vitro models of MM. M5 is proposed as a lead compound to be further optimized for improving its affinity toward MALAT1.

IRSOM2: a web server for predicting bifunctional RNAs.

Recent advances have shown that some biologically active non-coding RNAs (ncRNAs) are actually translated into polypeptides that have a physiological function as well. This paradigm shift requires adapted computational methods to predict this new class of 'bifunctional RNAs'. Previously, we developed IRSOM, an open-source algorithm to classify non-coding and coding RNAs. Here, we use the binary statistical model of IRSOM as a ternary classifier, called IRSOM2, to identify bifunctional RNAs as a rejection of the two other classes. We present its easy-to-use web interface, which allows users to perform predictions on large datasets of RNA sequences in a short time, to re-train the model with their own data, and to visualize and analyze the classification results thanks to the implementation of self-organizing maps (SOM). We also propose a new benchmark of experimentally validated RNAs that play both protein-coding and non-coding roles, in different organisms. Thus, IRSOM2 showed promising performance in detecting these bifunctional transcripts among ncRNAs of different types, such as circRNAs and lncRNAs (in particular those of shorter lengths). The web server is freely available on the EvryRNA platform: https://evryrna.ibisc.univ-evry.fr.

Investigating RNA-RNA interactions through computational and biophysical analysis.

Numerous viruses utilize essential long-range RNA-RNA genome interactions, specifically flaviviruses. Using Japanese encephalitis virus (JEV) as a model system, we computationally predicted and then biophysically validated and characterized its long-range RNA-RNA genomic interaction. Using multiple RNA computation assessment programs, we determine the primary RNA-RNA interacting site among JEV isolates and numerous related viruses. Following in vitro transcription of RNA, we provide, for the first time, characterization of an RNA-RNA interaction using size-exclusion chromatography coupled with multi-angle light scattering and analytical ultracentrifugation. Next, we demonstrate that the 5' and 3' terminal regions of JEV interact with nM affinity using microscale thermophoresis, and this affinity is significantly reduced when the conserved cyclization sequence is not present. Furthermore, we perform computational kinetic analyses validating the cyclization sequence as the primary driver of this RNA-RNA interaction. Finally, we examined the 3D structure of the interaction using small-angle X-ray scattering, revealing a flexible yet stable interaction. This pathway can be adapted and utilized to study various viral and human long-non-coding RNA-RNA interactions and determine their binding affinities, a critical pharmacological property of designing potential therapeutics.

Rolling circle extension-assisted loop-mediated isothermal amplification (Rol-LAMP) method for locus-specific and visible detection of RNA N6-methyladenosine.

N6-methyladenosine (m6A) is the most prevalent RNA modification in eukaryotic mRNAs. Currently available detection methods for locus-specific m6A marks rely on RT-qPCR, radioactive methods, or high-throughput sequencing. Here, we develop a non-qPCR, ultrasensitive, isothermal, and naked-eye visible method for m6A detection based on rolling circle amplification (RCA) and loop-mediated isothermal amplification (LAMP), named m6A-Rol-LAMP, to verify putative m6A sites in transcripts obtained from the high-throughput data. When padlock probes hybridize to the potential m6A sites on targets, they are converted to circular form by DNA ligase in the absence of m6A modification, while m6A modification hinders the sealing of padlock probes. Subsequently, Bst DNA polymerase-mediated RCA and LAMP allow the amplification of the circular padlock probe to achieve the locus-specific detection of m6A. Following optimization and validation, m6A-Rol-LAMP can ultra-sensitively and quantitatively determine the existence of m6A modification on a specific target site as low as 100 amol under isothermal conditions. Detections of m6A can be performed on rRNA, mRNA, lincRNA, lncRNA and pre-miRNA from biological samples with naked-eye observations after dye incubation. Together, we provide a powerful tool for locus-specific detection of m6A, which can simply, quickly, sensitively, specifically, and visually determine putative m6A modification on RNA.

LncReader: identification of dual functional long noncoding RNAs using a multi-head self-attention mechanism.

Long noncoding ribonucleic acids (RNAs; LncRNAs) endowed with both protein-coding and noncoding functions are referred to as 'dual functional lncRNAs'. Recently, dual functional lncRNAs have been intensively studied and identified as involved in various fundamental cellular processes. However, apart from time-consuming and cell-type-specific experiments, there is virtually no in silico method for predicting the identity of dual functional lncRNAs. Here, we developed a deep-learning model with a multi-head self-attention mechanism, LncReader, to identify dual functional lncRNAs. Our data demonstrated that LncReader showed multiple advantages compared to various classical machine learning methods using benchmark datasets from our previously reported cncRNAdb project. Moreover, to obtain independent in-house datasets for robust testing, mass spectrometry proteomics combined with RNA-seq and Ribo-seq were applied in four leukaemia cell lines, which further confirmed that LncReader achieved the best performance compared to other tools. Therefore, LncReader provides an accurate and practical tool that enables fast dual functional lncRNA identification.

ortho2align: a sensitive approach for searching for orthologues of novel lncRNAs.

BACKGROUND: Many novel long noncoding RNAs have been discovered in recent years due to advances in high-throughput sequencing experiments. Finding orthologues of these novel lncRNAs might facilitate clarification of their functional role in living organisms. However, lncRNAs exhibit low sequence conservation, so specific methods for enhancing the signal-to-noise ratio were developed. Nevertheless, current methods such as transcriptomes comparison approaches or searches for conserved secondary structures are not applicable to novel, previously unannotated lncRNAs by design. RESULTS: We present ortho2align-a versatile sensitive synteny-based lncRNA orthologue search tool with statistical assessment of sequence conservation. This tool allows control of the specificity of the search process and optional annotation of found orthologues. ortho2align shows similar performance in terms of sensitivity and resource usage as the state-of-the-art method for aligning orthologous lncRNAs but also enables scientists to predict unannotated orthologous sequences for lncRNAs in question. Using ortho2align, we predicted orthologues of three distinct classes of novel human lncRNAs in six Vertebrata species to estimate their degree of conservation. CONCLUSIONS: Being designed for the discovery of unannotated orthologues of novel lncRNAs in distant species, ortho2align is a versatile tool applicable to any genomic regions, especially weakly conserved ones. A small amount of input files makes ortho2align easy to use in orthology studies as a single tool or in bundle with other steps that researchers will consider sensible. ortho2align is available as an Anaconda package with its source code hosted at https://github.com/dmitrymyl/ortho2align .

A toolkit for the identification of NEAT1_2/paraspeckle modulators.

Paraspeckles are ribonucleoprotein granules assembled by NEAT1_2 lncRNA, an isoform of Nuclear Paraspeckle Assembly Transcript 1 (NEAT1). Dysregulation of NEAT1_2/paraspeckles has been linked to multiple human diseases making them an attractive drug target. However currently NEAT1_2/paraspeckle-focused translational research and drug discovery are hindered by a limited toolkit. To fill this gap, we developed and validated a set of tools for the identification of NEAT1_2 binders and modulators comprised of biochemical and cell-based assays. The NEAT1_2 triple helix stability element was utilized as the target in the biochemical assays, and the cellular assay ('ParaQuant') was based on high-content imaging of NEAT1_2 in fixed cells. As a proof of principle, these assays were used to screen a 1,200-compound FDA-approved drug library and a 170-compound kinase inhibitor library and to confirm the screening hits. The assays are simple to establish, use only commercially-available reagents and are scalable for higher throughput. In particular, ParaQuant is a cost-efficient assay suitable for any cells growing in adherent culture and amenable to multiplexing. Using ParaQuant, we identified dual PI3K/mTOR inhibitors as potent negative modulators of paraspeckles. The tools we describe herein should boost paraspeckle studies and help guide the search, validation and optimization of NEAT1_2/paraspeckle-targeted small molecules.

Regulation and roles of RNA modifications in aging-related diseases.

With the aging of the global population, accumulating interest is focused on manipulating the fundamental aging-related signaling pathways to delay the physiological aging process and eventually slow or prevent the appearance or severity of multiple aging-related diseases. Recently, emerging evidence has shown that RNA modifications, which were historically considered infrastructural features of cellular RNAs, are dynamically regulated across most of the RNA species in cells and thereby critically involved in major biological processes, including cellular senescence and aging. In this review, we summarize the current knowledge about RNA modifications and provide a catalog of RNA modifications on different RNA species, including mRNAs, miRNAs, lncRNA, tRNAs, and rRNAs. Most importantly, we focus on the regulation and roles of these RNA modifications in aging-related diseases, including neurodegenerative diseases, cardiovascular diseases, cataracts, osteoporosis, and fertility decline. This would be an important step toward a better understanding of fundamental aging mechanisms and thereby facilitating the development of novel diagnostics and therapeutics for aging-related diseases.

Distinct MUNC lncRNA structural domains regulate transcription of different promyogenic factors.

Many lncRNAs have been discovered using transcriptomic data; however, it is unclear what fraction of lncRNAs is functional and what structural properties affect their phenotype. MUNC lncRNA (also known as DRReRNA) acts as an enhancer RNA for the Myod1 gene in cis and stimulates the expression of other promyogenic genes in trans by recruiting the cohesin complex. Here, experimental probing of the RNA structure revealed that MUNC contains multiple structural domains not detected by prediction algorithms in the absence of experimental information. We show that these specific and structurally distinct domains are required for induction of promyogenic genes, for binding genomic sites and gene expression regulation, and for binding the cohesin complex. Myod1 induction and cohesin interaction comprise only a subset of MUNC phenotype. Our study reveals unexpectedly complex, structure-driven functions for the MUNC lncRNA and emphasizes the importance of experimentally determined structures for understanding structure-function relationships in lncRNAs.

Structural effects of m6A modification of the Xist A-repeat AUCG tetraloop and its recognition by YTHDC1.

The A-repeat region of the lncRNA Xist is critical for X inactivation and harbors several N6-methyladenosine (m6A) modifications. How the m6A modification affects the conformation of the conserved AUCG tetraloop hairpin of the A-repeats and how it can be recognized by the YTHDC1 reader protein is unknown. Here, we report the NMR solution structure of the (m6A)UCG hairpin, which reveals that the m6A base extends 5' stacking of the A-form helical stem, resembling the unmethylated AUCG tetraloop. A crystal structure of YTHDC1 bound to the (m6A)UCG tetraloop shows that the (m6A)UC nucleotides are recognized by the YTH domain of YTHDC1 in a single-stranded conformation. The m6A base inserts into the aromatic cage and the U and C bases interact with a flanking charged surface region, resembling the recognition of single-stranded m6A RNA ligands. Notably, NMR and fluorescence quenching experiments show that the binding requires local unfolding of the upper stem region of the (m6A)UCG hairpin. Our data show that m6A can be readily accommodated in hairpin loop regions, but recognition by YTH readers requires local unfolding of flanking stem regions. This suggests how m6A modifications may regulate lncRNA function by modulating RNA structure.

SMRT sequencing of the full-length transcriptome of Gekko gecko.

Tokay Gecko (Gekko gecko) is a rare and endangered medicinal animal in China. Its dry body has been used as an anti-asthmatic agent for two thousand years. To date, the genome and transcriptome of this species remain poorly understood. Here, we adopted single molecule real-time (SMRT) sequencing to obtain full-length transcriptome data and characterized the transcriptome structure. We identified 882,273 circular consensus (CCS) reads, including 746,317 full-length nonchimeric (FLNC) reads. The transcript cluster analysis revealed 212,964 consensus sequences, including 203,994 high-quality isoforms. In total, 111,372 of 117,888 transcripts were successfully annotated against eight databases (Nr, eggNOG, Swiss-Prot, GO, COG, KOG, Pfam and KEGG). Furthermore, 23,877 alternative splicing events, 169,128 simple sequence repeats (SSRs), 10,437 lncRNAs and 7,932 transcription factors were predicted across all transcripts. To our knowledge, this report is the first to document the G. gecko transcriptome using SMRT sequencing. The full-length transcript data might accelerate transcriptome research and lay the foundation for further research on G. gecko.

Temporal transcriptomic changes in long non-coding RNAs and messenger RNAs involved in the host immune and metabolic response during Toxoplasma gondii lytic cycle.

BACKGROUND: Long non-coding RNAs (lncRNAs) are important regulators of various biological and pathological processes, in particular the inflammatory response by modulating the transcriptional control of inflammatory genes. However, the role of lncRNAs in regulating the immune and inflammatory responses during infection with the protozoan parasite Toxoplasma gondii remains largely unknown. METHODS: We performed a longitudinal RNA sequencing analysis of human foreskin fibroblast (HFF) cells infected by T. gondii to identify differentially expressed long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs), and dysregulated pathways over the course of T. gondii lytic cycle. The transcriptome data were validated by qRT-PCR. RESULTS: RNA sequencing revealed significant transcriptional changes in the infected HFFs. A total of 697, 1234, 1499, 873, 1466, 561, 676 and 716 differentially expressed lncRNAs (DElncRNAs), and 636, 1266, 1843, 2303, 3022, 1757, 3088 and 2531 differentially expressed mRNAs (DEmRNAs) were identified at 1.5, 3, 6, 9, 12, 24, 36 and 48 h post-infection, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DElncRNAs and DEmRNAs revealed that T. gondii infection altered the expression of genes involved in the regulation of host immune response (e.g., cytokine-cytokine receptor interaction), receptor signaling (e.g., NOD-like receptor signaling pathway), disease (e.g., Alzheimer's disease), and metabolism (e.g., fatty acid degradation). CONCLUSIONS: These results provide novel information for further research on the role of lncRNAs in immune regulation of T. gondii infection.

exoRBase 2.0: an atlas of mRNA, lncRNA and circRNA in extracellular vesicles from human biofluids.

Extracellular vesicles (EVs) are small membranous vesicles that contain an abundant cargo of different RNA species with specialized functions and clinical implications. Here, we introduce an updated online database (http://www.exoRBase.org), exoRBase 2.0, which is a repository of EV long RNAs (termed exLRs) derived from RNA-seq data analyses of diverse human body fluids. In exoRBase 2.0, the number of exLRs has increased to 19 643 messenger RNAs (mRNAs), 15 645 long non-coding RNAs (lncRNAs) and 79 084 circular RNAs (circRNAs) obtained from ∼1000 human blood, urine, cerebrospinal fluid (CSF) and bile samples. Importantly, exoRBase 2.0 not only integrates and compares exLR expression profiles but also visualizes the pathway-level functional changes and the heterogeneity of origins of circulating EVs in the context of different physiological and pathological conditions. Our database provides an attractive platform for the identification of novel exLR signatures from human biofluids that will aid in the discovery of new circulating biomarkers to improve disease diagnosis and therapy.

Long Noncoding RNA MIR2187HG Suppresses TBK1-Mediated Antiviral Signaling by Deriving miR-2187-3p in Teleost Fish.

Long noncoding RNAs (lncRNAs) function as microregulatory factors that influence gene expression after a variety of pathogenic infections, and they have been extensively studied in the past few years. Although less attention has been paid to lncRNAs in lower vertebrates than in mammals, current studies reveals that lncRNAs play a vital role in fish stimulated by pathogens. Here, we discovered a new lncRNA, termed MIR2187HG, which can function as a precursor of a small RNA, miR-2187-3p, with regulatory functions in the miiuy croaker (Miichthys miiuy). Upon Siniperca chuatsi rhabdovirus (SCRV) virus infection, the expression levels of MIR2187HG were remarkably enhanced. Elevated MIR2187HG expression can act as a pivotally negative regulator that participates in the innate immune response of teleost fish to inhibit the intracellular TANK-binding kinase 1 (TBK1)-mediated antiviral signaling pathways, which can effectively avoid excessive immunity. In addition, we found that SCRV could also utilize MIR2187HG to enhance its own numbers. Our results not only provide evidence regarding the involvement of the lncRNAs in response to viruses in fish but also broaden our understanding of the function of lncRNAs as precursor microRNAs (miRNAs) in teleost fish for the first time. IMPORTANCE SCRV infection upregulates MIR2187HG levels, which in turn suppresses SCRV-triggered type I interferon production, thus promoting viral replication in miiuy croaker. Notably, MIR2187HG regulates the release of miR-2187-3p, and TBK1 is a target of miR-2187-3p. MIR2187HG could acquire from miR-2187-3p the function of inhibiting TBK1 expression and subsequently modulate TBK1-mediated NF-κB and IRF3 signaling. The collective results suggest that the novel regulation mechanism of TBK1-mediated antiviral response during RNA viral infection was regulated by MIR2187HG. Therefore, a new regulation mechanism for lncRNAs to regulate antiviral immune responses in fish is proposed.

Dynamic bulge nucleotides in the KSHV PAN ENE triple helix provide a unique binding platform for small molecule ligands.

Cellular and virus-coded long non-coding (lnc) RNAs support multiple roles related to biological and pathological processes. Several lncRNAs sequester their 3' termini to evade cellular degradation machinery, thereby supporting disease progression. An intramolecular triplex involving the lncRNA 3' terminus, the element for nuclear expression (ENE), stabilizes RNA transcripts and promotes persistent function. Therefore, such ENE triplexes, as presented here in Kaposi's sarcoma-associated herpesvirus (KSHV) polyadenylated nuclear (PAN) lncRNA, represent targets for therapeutic development. Towards identifying novel ligands targeting the PAN ENE triplex, we screened a library of immobilized small molecules and identified several triplex-binding chemotypes, the tightest of which exhibits micromolar binding affinity. Combined biophysical, biochemical, and computational strategies localized ligand binding to a platform created near a dinucleotide bulge at the base of the triplex. Crystal structures of apo (3.3 Å) and ligand-soaked (2.5 Å) ENE triplexes, which include a stabilizing basal duplex, indicate significant local structural rearrangements within this dinucleotide bulge. MD simulations and a modified nucleoside analog interference technique corroborate the role of the bulge and the base of the triplex in ligand binding. Together with recently discovered small molecules that reduce nuclear MALAT1 lncRNA levels by engaging its ENE triplex, our data supports the potential of targeting RNA triplexes with small molecules.

The lupus autoantigen La/Ssb is an Xist-binding protein involved in Xist folding and cloud formation.

Using the programmable RNA-sequence binding domain of the Pumilio protein, we FLAG-tagged Xist (inactivated X chromosome specific transcript) in live mouse cells. Affinity pulldown coupled to mass spectrometry was employed to identify a list of 138 candidate Xist-binding proteins, from which, Ssb (also known as the lupus autoantigen La) was validated as a protein functionally critical for X chromosome inactivation (XCI). Extensive XCI defects were detected in Ssb knockdown cells, including chromatin compaction, death of female mouse embryonic stem cells during in vitro differentiation and chromosome-wide monoallelic gene expression pattern. Live-cell imaging of Xist RNA reveals the defining XCI defect: Xist cloud formation. Ssb is a ubiquitous and versatile RNA-binding protein with RNA chaperone and RNA helicase activities. Functional dissection of Ssb shows that the RNA chaperone domain plays critical roles in XCI. In Ssb knockdown cells, Xist transcripts are unstable and misfolded. These results show that Ssb is critically involved in XCI, possibly as a protein regulating the in-cell structure of Xist.